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Abstract 

Regulation of mRNA decay is a critical component of global cellu- 
lar adaptation to changing environments. The corresponding changes in 
mRNA lifetimes can be coordinated with changes in mRNA transcrip- 
tion rates to fine-tune gene expression. Current approaches for measuring 
mRNA lifetimes can give rise to secondary effects due to transcription in- 
hibition and require separate experiments to estimate changes in mRNA 
transcription rates. Here, we propose an approach for simultaneous deter- 
mination of changes in mRNA transcription rate and lifetime using regula- 
tory small RNAs to control mRNA decay. We analyze a stochastic model 
for coupled degradation of mRNAs and sRNAs and derive exact results 
connecting RNA lifetimes and transcription rates to mean abundances. 
The results obtained show how steady-state measurements of RNA levels 
can be used to analyze factors and processes regulating changes in mRNA 
transcription and decay. 

Cellular adaptation to changing conditions is critically dependent on pro- 
cesses that enable rapid responses to environmental fluctuations. While consid- 
erable research has focused on changes in transcription, research over the past 
several years has demonstrated that control of mRNA decay plays an increas- 
ingly important role in cellular responses [TJ [5J |3J. Correspondingly, there is 
significant interest in understanding factors and processes which govern regula- 
tion of mRNA decay. 

The traditional approach for measuring mRNA lifetimes involves quantifica- 
tion of mRNAs remaining at different times following inhibition of transcription, 
e.g., by the addition of rifampicin [2J. This procedure requires multiple mea- 
surements during time intervals of the order of the mRNA lifetime, hence high 
temporal resolution is required for short-lived mRNAs. More significantly, the 
procedure for inhibition of transcription can give rise to secondary effects which 
influence mRNA decay j2j, hence it is of interest to consider alternative ap- 
proaches. In the following, we outline a novel proposal for quantifying mRNA 
decay. 
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Figure 1: Panel A: Kinetic scheme for mRNA and sRNA production and decay 
including coupled mutual degradation. Panel B: The proposed setup involves 
steady-state measurements for three strains: A(sRNA), WT and A (mRNA). 
Unregulated steady-state mean levels of mRNAs and sRNAs along with regu- 
lated levels of these molecules are measured. The measured quantities allow 
determination of the average mRNA transcription rate k m and decay rate T m 
relative to the sRNA production rate k s . If k s is held fixed, and the conditions 
are varied, the proposed scheme leads to simultaneous determination of fold- 
changes in the rate of transcription and the rate of mRNA decay. Note that the 
mRNA/sRNA interaction parameter can be arbitrary. 



A naturally occurring process which regulates turnover for many bacterial 
mRNAs involves interactions with regulatory small RNAs (sRNAs). Experi- 
ments have shown that post-transcriptional regulation of gene expression can 
occur via binding and subsequent coupled degradation of the mRNA and regu- 
latory sRN A which occurs rapidly |10[ I17j . The corresponding coarse-grained 
model of gene expression has been analyzed by several groups QTJ Q2J [13] 
and shown to be consistent with experimental observations [S]. We will show 
that the kinetic scheme for this mode of regulation leads to an exact mathemat- 
ical result relating RNA decay rates. The obtained result is valid even in the 
presence of large fluctuations that are typical for low abundance mRNAs. 

Our analysis considers a generalized stochastic model for regulation by sR- 
NAs as illustrated in Fig. 1 . The complexity of processes leading to transcription 
suggests that RNA synthesis in many instances is not adequately modeled as 
a Poisson process pj Q3]. Hence, we model mRNA and sRNA production by 
arbitrary stochastic processes with mean arrival rates k m and k s respectively. 
Degradation of a mRNA and a sRNA can occur either independently (with con- 
stant probability per unit time — and — respectively) or through the coupled 
degradation process. 

The experimental setup in our proposal for quantifying mRNA decay is as 
follows. Consider three different strains as shown in Fig. IB: two unregulated 
strains (i.e., with cither sRNA or mRNA deleted) and the wild type (WT) strain. 
In the WT strain, both mRNA and sRNA are present and regulate each other. 
In steady state, we derive the following exact relations (Appendix) connecting 
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mRNA/sRNA lifetimes and transcription rates to the mean abundances 
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Here (to) and (s) are the mean mRNA(sRNA) abundances in strains lacking 
the sRNA(mRNA), and (to) and (s) are the mean mRNA and sRNA levels in 
WT strain where both are present. 

The above relations suggest an alternative approach (Fig. IB) for quantify- 
ing decay times for mRNAs that either have a naturally occurring small RNA 
regulator or for which an antisense RNA regulator can be designed. Consider 
an experimental setup expressing the sRNA from an inducible promoter such 
that its transcription rate is primarily controlled by inducer concentration. The 
mean transcription rate k s can, in principle, be determined using single-molecule 
methods [18]. Now the basic parameters for the coupled system are k m , k s , r m 
and t s . If k s is known, then the values of the other parameters can be deter- 
mined using experimental measurements of (to), (s), (to) and (s) in combination 
with equations Eqns. 1-3 above. Alternatively, experiments can be designed to 
keep k s fixed while factors regulating mRNA decay are changed e.g., by dele- 
tion of a protein known to play a role in mRNA decay. The above equations 
can be used to simultaneously determine the corresponding fold-changes of the 
mRNA/sRNA lifetimes and the mean mRNA transcription rate. 

The proposed approach can be used to address several important questions 
of current interest, some of which are highlighted in the following. By targeted 
mutagenesis of specific mRNA sequence elements, the induced fold-change in 
mRNA lifetime, as well as the corresponding change in the transcription rate k m , 
can be determined using the same experimental setup. This is an important fea- 
ture, given that recent experiments have observed coordination between changes 
in transcription and changes in mRNA degradation |16j . Quantifying the change 
in mRNA lifetimes induced by mutations to different components of cellular 
degradation pathways can address such issues as the role of polyadenylation in 
mRNA decay Hj. It would also be of interest to design high-throughput ex- 
periments for different mRNAs which are regulated by corresponding antisense 
RNAs, all of which are expressed from identical inducible promoters and thus 
have the same k s . The proposed procedure can then be used for genome- wide 
determination of relative transcription rates and lifetimes of mRNAs. These 
effective parameters, in turn, serve as critical inputs to systems-level models of 
cellular processes [T3] , 

In summary, we have proved an exact relation for a nonlinear stochastic 
model of cellular post-transcriptional regulation. The derived results suggest 
a novel procedure for simultaneous determination of mRNA production rates 
and mRNA lifetimes. While the focus was on bacterial mRNAs, the procedure 
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outlined can also be applied to higher organisms and used to systematically ex- 
plore the sequence determinants and processes involved in regulation of mRNA 
decay. 

Appendix 

The proposed experimental setup involves measurements of mRNA/sRNA abun- 
dances for three strains, the sRNA deleted A(sRNA) strain, mRNA deleted 
A (mRNA) strain, and the wild type (WT) strain (both species are present 
and regulate each other.) The regulation model (with constant creation rates 
for mRNA and sRNA) has been analyzed by several groups [5] [TTJ |T3] and is 
summarized in the reaction scheme described in Fig.l. 

Recent experiments which provide evidence for transcriptional bursting |14j 
point towards the need to go beyond the Poisson process as a model for RNA 
synthesis. Accordingly, we take both mRNA and sRNA creation events as 
arbitrary stochastic processes. Degradation of RNA is assumed to be a Poisson 
process with rate t x x , where t x is the mean lifetime of irRNA , x — {to, s}. In 
the WT strain, mRNA and sRNA undergo a coupled degradation process. We 
assume that this process is symmetric with respect to the number of mRNA 
and sRNA molecules involved, e.g., M + S — > 0. 

Let us choose a particular realization of the system evolution during time 
interval t = [0, T]. For large values of T, we derive 

x(T) - x(Q) = C X (T) - Y(T) - r- 1 [ dtx(t), (3) 

Jo 

where x(t) is the number of molecules of the species x — {to, s} at the time t. 

In Eq.[3] C x (t) is the total number of molecules of the species x created during 
system evolution until time T, and Y(T) is the total number of molecules of 
either species that is mutually degraded within the time interval [0, T]. Finally, 
using the law of large numbers, the number of molecules degraded naturally in 
[0, T] is given by the last term in the Eq.[3] 

Dividing both sides of Eq.[3]by T and taking a limit T — > oo we obtain 

r x{T)-x(0) Y(T) 

hm ^ — — — k x - t„ (x) - hm , (4) 

where (x) is average number of molecules in the system. We also defined k x as 
a mean arrival rate of the species x — {to, s}. The limit on the left hand side of 
Eq. [4] vanishes in the case of finite degradation rates r" 1 (number of molecules 
at any time is finite.) Note that Y(T) is an extensive quantity (it is monotonic 
increasing function of T) and therefore, the limit on the right hand side of Eq.|4] 
is finite. 
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Hence, we derive 



k m -Tj-{m) - lim = 0, 

1 — *OQ 1 

fc,-r i - 1 (S)- lim ^P- = 0, (5) 
which immediately yields the following expression 

fc rn - T" 1 (m) = fc s - T7 1 (s) . (6) 

In the unregulated case Y(T) = for any T, since one of the RNA species 
is deleted and there is no coupled degradation. In this situation one gets 

= fc m -r- 1 (TO),0 = fc s -r s - 1 (s), (7) 

where (x), x — {to, s} are the average number of molecules during unregulated 
system evolution. Combining the set of equations above with Eq.[6] we derive 
the results in the main text Eqns 1-3. We note that the derived results are 
valid even if the binding of mRNA and sRNA is taken to be reversible and the 
lifetime of the mRNA-sRNA complex is finite. Finally, the time average can be 
replaced by the ensemble average in the steady state. 

We have validated the derived results using stochastic simulations based 
on the Gillespie algorithm 0]. Production of RNA molecules was taken to 
occur in transcriptional bursts [2] i.e. each burst corresponds to the arrival of 
a random number of RNAs drawn from a geometric distribution, conditional 
on the production of at least 1 RNA molecule [5 . The waiting-time between 
bursts was a random variable drawn from exponential or Gamma distributions. 
As expected, the results from the simulations were in excellent agreement with 
the derived analytical results. 
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